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Abstract 

Solutions are constructed to the quantum constraints for planar 
gravity (fields dependent on z and t only) in the Ashtekar complex 
connection formalism. A number of operators are constructed and 
applied to the solutions. These include the familiar ADM energy and 
area operators, as well as new operators sensitive to directionality 
(z+ct vs. z-ct dependence). The directionality operators are quan- 
tum analogs of the classical constraints proposed for unidirectional 
plane waves by Bondi, Pirani, and Robinson (BPR). It is argued 
that the quantum BPR constraints will predict unidirectionality re- 
liably only for solutions which are semiclassical in a certain sense. 
The ADM energy and area operators are likely to have imaginary 
eigenvalues, unless one either shifts to a real connection, or allows 
the connection to occur other than in a holonomy. In classical the- 
ory, the area can evolve to zero. A quantum mechanical mechanism 
is proposed which would prevent this collapse. 

PACS categories: 04.60, 04.30 



I Introduction: Classical Radiation Criteria 

The connection-triad variables introduced by Ashtekar have sim- 
plified the constraint equations of quantum gravity; further, these 
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variables suggest that in the future we may be able to reformulate 
gravity in terms of non-local holonomies rather than local field op- 
erators 0, [| f3j. However, the new variables are unfamiliar, and it 
is not always clear what they mean physically and geometrically. In 
particular, it is not clear what operators or structures correspond 
to gravity waves. Although the quantum constraint equations are 
much simpler in the new variables, and solutions to these equations 
have been found P, |5| f|, it is not clear whether any of these solu- 
tions contain gravitational radiation. 

This is the fifth of a series of papers which search for operator 
signatures for gravitational radiation by applying the Ashtekar for- 
malism to the problem of plane gravitational waves. Paper I in 
the series |7j constructed classical constants of the motion for the 
plane wave case, using the more familiar geometrodynamics rather 
than Ashtekar connection dynamics. Papers II and III switched to 
connection dynamics and proposed solutions to the quantum con- 
straints || |9|] . The constraints annihilate the solutions of II except 
at boundary points, and annihilate the solutions of III everywhere. 
Paper IV constructs an operator Lz which measures total intrinsic 



spin around the z axis |10[. The present paper proposes operator 
signatures which are sensitive to the directionality of gravitational 
radiation (z -ct vs. z+ct dependence) and applies those operators 
(as well as the spin, energy, and area operators) to the solutions 
constructed in II and III. 

It is not easy to detect the presence of radiation, even when the 
problem is formulated classically, using the more familiar metric 
variables. One would like to define gravitational radiation using an 
energy criterion (as: radiation is a means of transporting energy 
through empty space ...). Gravitational energy is notoriously diffi- 
cult to define, however, since there is no first-order-in-derivatives-of 
-the -metric quantity which is a tensor. Accordingly, in the period 
1960-1970 several authors developed an algebraic criterion involving 
transverse components of a second-order quantity, the Weyl tensor 
11, 12, 13, H] . To use the criterion, one needs to know which 



directions are "transverse"; hence the criterion is most useful when 
the direction of propagation is clear from the symmetry: e.g. ra- 
dial propagation (for spherical symmetry) or z-axis propagation (for 
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planar symmetry, the case studied in the present paper). The pla- 



nar metrics considered here [HR Il6l admit two null vectors k and 1 



which have the right hypersurface orthogonality properties to be the 
propagation vectors for right-moving (k) and left-moving (1) grav- 
itational waves along the z axis, so that the propagation direction 
is especially easy to identify. The Weyl criterion is derived and 
discussed in Appendix D. 

A second, more group-theoretical criterion was developed by Bondi, 
Pirani, and Robinson (BPR) JL7[]. It is applicable when the plane 
wave is unidirectional, that is, when the wave is either right-moving 
(depending only on z-ct) or left-moving (depending only on z+ct). 
The unidirectional case is especially intruguing. It is relatively 



simple, since no scattering occurs [[18], [L9|]. Nevertheless the full 
complexity of gravity is already present; the unidirectional case is 
not simply waves propagating in an inert background. In particu- 
lar, no one has been able to cast the Hamiltonian into a free-field 
form in terms of variables (7^, q J ) which commute in the canoni- 
cal [7Ti, qj] = — iM- manner characteristic of non-interacting, non- 
gravitational theories. Also, the BPR criterion for unidirectional 
radiation requires three amplitudes to vanish, rather than the two 
one would naively expect from counting the two polarizations associ- 
ated with unidirectional radiation. I shall argue that the remaining 
vanishing amplitude represents a constraint on the background ge- 
ometry, a constraint which must be satisfied in order for the waves 
to propagate without backscattering off the background. The BPR 
group is derived and discussed in section II. 

Note that one criterion, that based on the Weyl tensor, is rele- 
gated to an appendix; while the BPR criterion is discussed in the 
body of the paper. I have done this primarily because the BPR 
amplitudes are much simpler than the Weyl amplitudes. It is not 
possible to ingnore the Weyl amplitudes completely, however; they 
are central to the literature of the 60's. Further, expressions which 
appear complex at one time may appear simple at a later time. At 
one time the traditional scalar constraint was thought to be too 
complex because it contained a factor of l/^ s /g■ Then Thiemann 
proposed a regularization of this constraint which actually requires 
that factor [EQ]. Similarly, the present "complexity" of the Weyl 
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amplitudes may disappear once a quantum regularization is con- 
structed. 

It is of some interest to reexpresss the two classical criteria, Weyl 
and BPR, in the Ashtekar language, even if one does not go on 
to consider the quantum case. However, one would really like to 
construct from each classical criterion a corresponding quantum op- 
erator. I construct such operators in section III and argue that these 
operators are reliable only when acting on wavefunctionals for states 
which are semiclassical, in a sense to be defined in section III. 

In section IV, I construct additional solutions to the quantum 
constraints. In section V, I apply the BPR quantum operators to 
the wavefunctional solutions obtained in papers II-III, as well as 
the new solutions constructed in section IV. Also in section V, I 
apply the ADM energy operator to the solutions, as well as the area 
operator and the operator L z for total intrinsic spin. 

There are three appendices. Two of the appendices (A and C) 
cover calculational details and the details of the Weyl criterion not 
used in the body of the paper. Appendix B considers the ADM 
energy. There is a modest surprise here: normally the ADM energy 
is considered to be given by the surface term in the Hamiltonian; but 
in the quantum case it is possible for the volume term to contribute 
also. 

My notation is typical of papers based upon the Hamiltonian 
approach with concomitant 3 + 1 splitup. Upper case indices A, 
B, . . .,1, J, K, ... denote local Lorentz indices ("internal" SU(2) 
indices) ranging over X, Y, Z only. Lower case indices a, b, . . ., i, j, 
. . . are also three- dimensional and denote global coordinates on the 
three-manifold. Occasionally the formula will contain a field with a 
superscript (4), in which case the local Lorentz indices range over 
X, Y, Z, T and the global indices are similarly four-dimensional; or 
a (2), in which case the local indices range over X, Y (and global 
indices over x, y) only. The (2) and (4) are also used in conjunction 
with determinants; e. g., g is the usual 3x3 spatial determinant, 
while ( 2 )e denotes the determinant of the 2x2 X, Y subblock of the 
triad matrix ■ I use Levi- Civita symbols of various dimensions: 
^txyz = cxyz = e XY = +1- The basic variables of the Ashtekar 
approach are an inverse densitized triad and a complex SU(2) 
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connection A A . 

E\ = eel; (1) 
[El, A*] = M(x-x')5X (2) 
The planar symmetry (two spacelike commuting Killing vectors, 



d x and d y in appropriate coordinates) allows Husain and Smolin [^1 
to solve and eliminate four constraints (the x and y vector constraint 
and the X and Y Gauss constraint) and correspondingly eliminate 
four pairs of (E^, A^) components. The 3x3 matrix then assumes 
a block diagonal form, with one lxl subblock occupied by E\ plus 
one 2x2 subblock which contains all the "transverse" E^, that is, 
those with a = x,y and A = X,Y. The 3x3 matrix of connections A^ 
assumes a similar block diagonal form. None of the surviving fields 
depends on x or y. 

The local Lorentz indices are vector rather than spinor; strictly 
speaking the internal symmetry is 0(3) rather than SU(2), gauge- 
fixed to 0(2) rather than U(l). Often it is convenient to shift to 
transverse fields which are eigenstates of the surviving gauge invari- 
ance 0(2): 

E% = [E|±iE^]/V2, (3) 

where a = x,y; and similarly for A^. 

In papers I-III I use the letter H to denote a constraint (scalar, 
vector, or Gauss). In the present paper I adopt what is becoming a 
more common convention in the literature and use the letter C to 
denote a constraint, while reserving the letter H for the Hamiltonian. 
The quantity denoted Cs in the present paper is identical to the 
constraint denoted Hg in papers II-III. This convention underscores 
the fact that every gravitational theory has constraints, but not 
every gravitational theory has a Hamiltonian. 

In three spatial dimensions it is usual to place the boundary 
surface at spatial infinity. Bringing the surface at infinity in to 
finite points is a major change, because at infinity the metric goes 
over to flat space, and flat space is a considerable simplification. In 
the present case (effectively one dimensional because of the planar 
symmetry) the space does not become flat at z goes to infinity, and 
nothing is lost by considering an arbitrary location for the boundary 
surface. The "surface" in one dimension is of course just two points 
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(the two endpoints of a segment of the z axis). The notation Zb 
denotes either the left or right boundary point zi or z r , z\ < z < z r . 
The result that the space does not become flat as z goes to infinity 
was established in paper II. Note that this result agrees with one's 
intuition from Newtonian gravity, where the potential in one spatial 
dimension due to a bounded source does not fall off, but grows as z 
at large z. 

If a certain solution does not satisfy the Gauss constraint (or 
other constraint) at the boundary, this does not mean that neces- 
sarily there is something wrong with the solution. In classical theory 
the solutions satisfy the constraints everywhere. In quantum the- 
ory, however, when the constraints are imposed after quantization, 
in the Dirac manner, it is only necessary that the smeared constraint 
annihilate the solution: 

/ dz5N(z)C(z)^ = 0. (4) 

The expression 5NC generates a small change SN in the Lagrange 
multiplier N. If N obeys a boundary condition of the form N— > 
constant at boundaries z b , then eq. (Q) must respect this boundary 
condition, which means 

5N(z b ) = 0. (5) 

Eq. (H) implies that C(z h )ip does not have to vanish. A statement 
that "this solution does not obey the constraint at the boundaries" 
does not mean necessarily that the solution is flawed. 

II Bondi-Pirani-Robinson Symmetry 

Bondi, Pirani, and Robinson (BPR) argue that the metric of a unidi- 
rectional plane gravitational wave should be invariant under a five- 
parameter group of symmetries. Their argument proceeds essen- 
tially as follows. First they point out that a plane electromagnetic 
wave moving in the +z direction is invariant under a five parameter 
group. (Besides the obvious <9 X , d y , and d v symmetries, there are 
two "null rotations" which rotate the v = (t + z)/v^2 direction into 
the x or y direction.) Then for gravitational plane waves they con- 
struct five Killing vectors which have the same Lie algebra as the 
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corresponding Killing vectors for the electromagnetic case. (More 
precisely they construct ten Killing vectors, one set of five for ct+z 
waves, and a similar set of five for ct-z waves.) 

This section constructs the five ct-z vectors, then imposes the 
usual symmetry requirement that the Lie derivative of the basic 
Ashtekar fields must vanish in the direction of the Killing vectors. 
In this way one finds that the fields must obey certain constraints; 
the section closes with a discussion of the physical meaning of these 
constraints in the classical theory. 

It is convenient to do the proofs in a gauge which has been sim- 
plified as much as possible using the d v symmetry, then afterwards 
transform the results to a general gauge. The plane wave metrics 
we consider here possess two hypersurface orthogonal null vectors; 
if the two hypersurfaces are labeled u = constant and v = constant 
(u and v = (ct ± z)/ y/2), then one can always transform the metric 
to a conformally flat form in the (z,t) sector by using u and v as 
coordinates fl6fl : 

(2) 

ds 2 = -2dudvf(u, v) + gabdx a dx b . (6) 

The sums over a,b,c,. . . extend over x,y only. If one now invokes 
the symmetry under d v , then f(u,v) depends on u only, and one can 
remove the function f by transforming to a new u coordinate. In 
this gauge, Rosen gauge fl22| , the metric is (not just conformally 
flat, but) flat in the (z,t) sector and non-trivial only in the (x,y) 
sector. 

(2) 

ds 2 = -2dudv + £ g ab dx a dx b . (7) 
In Rosen gauge, the Killing vectors are d x , d y , d v , and 

The constraints imposed by the first three <9 X , <9 y , and d v Killing 
vectors are satisfied already, because of the choice of gauge. I now 
work out the constraints which the last two Killing vectors, eq. (||), 
impose on the Ashtekar variables (in Rosen gauge first, than in a 
general gauge). I summarize the highlights of the calculation in this 
section, and move the algebraic details to appendix A. 
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It is necessary to calculate the symmetry constraints on the four- 
dimensional tetrads and Ashtekar connections first, since the Killing 
vectors are intrinsically four-dimensional, then carry out a 3+1 de- 
composition to obtain the constraints on the usual three dimensional 
densitized triad and connection. At the four- dimensional level, the 
three local Lorentz boosts have been gauge-fixed by demanding that 
three of the tetrads vanish: 

e 4 M = 0, M = space. (9) 

The gauge condition of eq. (Bh is the standard choice, used with 



all metrics p3[ . In addition, for the special case of the plane wave 
metric, the gauge-fixing of the XY Gauss constraint and xy spatial 
diffeomorphism constraints imply that four more tetrads vanish |]2l] . 



e z x = e z Y = el = el = 0. (10) 

At the four- dimensional level, the requirement of vanishing Lie 
derivative in the direction of the Killing vector gives 

= £ A 0Ae?-fyre?-L5e?; (n) 

= £\9fAL J + d a e (4) A< J 

+C l v ^Ai J + C J y ^k l l' - d Q C u . (12) 

These equations are not quite the usual Lie derivatives because of 
the the L and C terms. L and C are local Lorentz transformations. 
If £ = <9 X , d y , or d v , no L or C is required. If £ = one of the two 
Killing vectors defined at eq. (§), then a Lorentz transformation L 
is required in eq. (0); otherwise the symmetry destroys the gauge 
conditions, eqs. @ and (0). Then for consistency, ^A in eq. (]12f) 
must undergo the same Lorentz transformation; since ^A is self- 
dual, the Lorentz transformation C in eq. ([12]) must be the self-dual 
version of the Lorentz transformation L: 

2£ IJ = L IJ + i5(e IJ MN /2e TXY z)L MN . (13) 

The phase £/ e TXYZ = ±1 is the duality eigenvalue which determines 
whether the theory is self-dual or anti-self-dual. Because I include 
the extra factor of exxYZ, eq. (|13|) contains two factors of e, so is 
independent of one's choice of phase for this quantity. After the 
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four-dimensional theory is rewritten in 3+1 form, all results will 
depend only on 5. In the body of the paper I choose 5 = +1, but 
appendix A indicates what happens for the opposite choice 5 = — 1. 

It is a straightforward matter to determine the Lorentz transfor- 
mation L which will preserve the gauge conditions of eqs. (§) and 
(|TT]|), then solve eqs. (|TT|) and flT2|). This is done in appendix A. From 
eq. (|A23|) in that appendix, eqs. (|TTD and (|12|) imply the following 



constraints on the connection A. 

= AJ; 

= -A+ + 2"Re"A+, 

(right-moving) (14) 

where a = x,y only. The connection A is now the usual 3+1 con- 
nection, not the four- dimensional connection ^ 4 ^A. Re A^ without 
the quotes is the usual real part, containing no factors of i, while 

"Re"A+ = (ReA* + iReAj)/V2. (15) 

"Re" A contains a factor of i, because of the i in the definition of 
the (X ± iY) 0(2) eigenstates, and is no longer real. If one writes 
out the "Re" and "Im" parts of A^ it is easy to see that the two 
constraints of eq. (|14]) are just the complex conjugates of each other. 
To obtain the constraints for left-moving waves, interchange + and 
- in eq. ([Lj). 

Eq. (|I|]) can be interpreted physically by using the classical equa- 
tions of motion to prove theorems about the spin behavior of the 
BPR fields. Again, the required calculations are done in an appendix 
(appendix B); and this section summarizes the main conclusions. 

To interpret the spin content of the four amplitudes which van- 
ish, it is better to work with the four combinations A~ and Eg 
[-A+ + 2 "Re" A + ]. (One can always recover the original four ampli- 
tudes from these four, because one can always invert the 2x2 matrix 
formed from the transverse components of the densitized triad, E^ 



with B = X,Y and b = x,y.) Eq. (|B4j ) of appendix B expresses the 



total spin angular momentum Lz of the gravitational wave in terms 
of these combinations. 

L z = i|dzKe +x E a _[A- + (A--2ReA-) 



+e y _e_ x E a + [A+ + (A+ - 2ReA+)} - (x<-y), (16) 

where eA a and e^ are triad and inverse triad fields, respectively. Out 
of the four possible combinations A~ and Eg [- A+ + 2 "Re" A+] , 
only two combinations E 3 ^ A~ and E^ [-A+ + 2 "Re" A+] contribute 
to the spin angular momentum. Since these are the two amplitudes 
with 0(2)_ helicity ± 2 in the local Lorentz frame, it is natural to 
interpret Efj.A^ as an amplitude for a wave having helicity ±2. Both 
helicity ±2 combinations must vanish in order to eliminate the two 
polarizations moving in the ct + z direction. 

Two more, helicity zero combinations, E^A^ and EL[— A+ + 



2"Re"A+], also contain the fields of eq. ([14]). How does one interpret 
these helicity zero amplitudes? 

Using the Gauss constraint plus the classical equations of motion 
in a conformally flat gauge, one can prove that these two constraints 
collapse to a single constraint, eq. ( |B10| ) of Appendix B. 

= (d t + d z )E z z . (17) 

El = e|e = We, where ^e is the determinant of the 2x2 transverse 
sector of the triad matrix, a scalar function of (z,t). This function 
would seem to characterize the background geometry, rather than 
the wave. In general relativity, however, "background" and wave are 
inseparable, in the sense that the "background" is not inert. The 
wave will scatter off the background, in general, unless it obeys the 
constraint given by eq. ([17]) • 

This completes the survey of the constraints predicted by BPR 
symmetry, and their physical interpretation in the classical context. 
In the next section these classical expressions are promoted to quan- 
tum operators. 



Ill The Transition from Classical to Quantum 
Criterion 

This section lists four issues which arise when converting a classical 
expression into a quantum criterion. I summarize and discuss each 
issue, then show the application to the BPR criteria. 
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A Factor ordering and regularization 



I choose a factor ordering which is natural and simple within the 
complex connection formalism. The ordering is "functional deriva- 



tives to the right" p4j . That is, I quantize in a standard manner, by 
replacing one half the fields by functional derivatives, 

E| _> -M/5A Z Z ; 

(for a = x,y and A = X,Y), (18) 

and then order the functional derivatives to the right in every oper- 
ator or constraint. Regularization, if needed, is via point splitting 



21]. This approach has the virtue of consistency, since I have used 



it in two previous papers on quantization of plane waves. 



B Semiclassicality 

An example from QED will be helpful in explaining what is meant 
by semiclassicality. The BPR criteria are essentially field strengths 
for waves moving in a given direction with a given polarization. An 
analogous quantity from flat space QED is 

T = J^m,X (19) 

Here is the self-dual QED field strength and (k,l,m,m ) is the 
usual flat space null tetrad: k and 1 are null vectors with space 
components along ± z ; m and m are transverse polarization vectors 
along (x ± iy) / y/2. Classically, the criterion for absence of radiation 
along k with polarization m is T = 0. The corresponding quantum 
criterion for absence of radiation, obtained after replacing classical 
A fields by quantum operators A is not 

Pip = 0, (wrong) (20) 

but rather 

(semicl] JF|semicl) = 0. (21) 

Since the QED field strength T contains creation as well as annihi- 
lation operators, it cannot annihilate any state, and eq. ( p0|) is too 
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strong. The classical statement T = merely implies the existence 
of a corresponding semiclassical state such that eq. ( pT| ) holds. I have 
deliberately used the term "semiclassical" rather than "coherent' to 
describe the state in eq. fl2"T|), because the latter term conventionally 
denotes a state which is an eigenfunction of the annihilation opera- 
tor, and annihilation operators usually are not available in quantum 
gravity. While annihilation operators may not exist, certainly semi- 
classical states will, because of the correspondence principle. 

Clearly the criterion eq. (|2lD is more difficult to apply than 
eq. (pUD, but let us survey the damage; the situation may not be 
hopeless. To define "semiclassical" without invoking coherent states 
or annihlation operators, one can study 

F|semicl) = f(z)|semicl) + [remainder), (22) 

where F stands for a typical BPR field. |semicl) is normed to unity, 
although | remainder) need not be. I define semiclassical, not by re- 
quiring f(z) to be large, but rather by requiring the [remainder) to 
be small. If I require f(z) to be large (perhaps reasoning that "clas- 
sical" means large quantum numbers) then I exclude the vanishing 
amplitude case, f(z) = 0, where |semicl) is the vacuum with respect 
to a given radiation mode. The vacuum is a well-defined state, clas- 
sically, and one expects it to have a quantum analog. For |semicl) to 
approximate a classical state, therefore, it is not necessary that f(z) 
be large; only that the fluctuations away from this state be small. 
These fluctuations are measured by 

(F f F) - (Ff)(F) = (remainder | remainder) < (l p /l) 2 /l 2 . (23) 

Here the fluctuations are assumed to be small compared to the size 
of typical matrix elements (f |F|i) « l p /l 2 one gets when |i) and (f | 
are few-graviton states and F is a canonical degree of freedom in 
the linearized theory. l p is the Planck length, and 1 is a typical 
length or wavelength. (To modify this discussion for the Weyl fields 
of appendix D, which are classical dimension l/l 2 , replace l p /l 2 by 
l p /l 3 . From now on ( ) is understood to indicate an average over a 
semiclassical state, unless explicitly indicated otherwise.) 

One might ask what is meant by a typical length 1, in a quantized 
and diffeomorphism invariant theory where no background metric is 
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available. In such a theory, even in the absence of a background met- 



ric, length, area, and volume operators can be defined pR E6|. 120 



and the eigenvalues of these geometric operators are dimensionless 
functions of spins ji, times factors of l p to give the correct dimension. 
The spins label the irreducible representations of SU(2) associated 
with each holonomy (if the wavefunctional is in connection repre- 
sentation, a product of holonomies) or associated with each edge 
of a spin network (if the wavefunctional is a spin network state). 
Thus one expects 1 = l (ji)l p , lo dimensionless and 3> 1. (In the 
planar case, SU(2) is gauge-fixed to 0(2) and presumably the SU(2) 
eigenvalues j will be replaced by the 0(2) eigenvalues m = spin an- 
gular momentum component along z.) Evidently, then, to check 
that eq. ( p3|) is satisfied, one should apply the geometric operators 
to the state first, in order to estimate 1. One also needs the measure 
in Hilbert space; but in favorable situations one might be able to 
tell that [remainder) is small simply by inspection of eq. (p2"|) . 

Since the criterion of eq. (|23| ) is an inequality, it cannot be used 
to draw sharp distinctions between states. For example, in the lin- 
earized limit, if |N) denotes an eigenstate of the number operator 
having N quanta of a given polarization and direction, then f(z) will 
be zero, while the norm eq. (p3|) will be order N (l p /l) 2 /l 2 . There 
will be uncertainties in estimating 1, so that the criterion cannot dis- 
tinguish sharply between the vacuum state and a number eigenstate 
having small occupation number N. 



C Non-polynomiality 

The BPR operators, eq. (0), occur in complex conjugate pairs, and 
one member of the BPR pair involves Re A^, which is a known, 
but non-polynomial function of E^. In particular, Re contains 
factors of l/^ 2 -E , where ^-E is the 2x2 determinant formed from 
the E^ with internal indices A = X,Y and global indices a = x,y. I 
have dealt with a similar operator, 1/E|, in a previous paper; but 
would just as soon not do so here. 

One can use the fact that the BPR constraints come in complex 
co jugate pairs, plus semiclassicality, to prove the following theo- 
rem (and then one uses the theorem to avoid dealing with the non- 
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polynomiality). Theorem: 



(-A+ + 2ReA+> = (A a >, (24) 

and similarly for the other BPR pair. This result is just what one 
would expect from the corresponding result for expectation values 
of complex operators in ordinary quantum mechanics (for instance 
(p + iq)* — (p — iq)) except that here the basic operators are not 
Hermitean, so that the proof is slightly longer. Proof: Expand out 
the operators using 

A± = (A?±iA*)/V2 

= U(5/5E^±i5/6E^)/V2. (25) 

Integrate by parts the functional derivative on the left side of eq. (|24"D , 
using 

^ * h8ip/8E% = J /i[(-h)<V£ElVHV> + J flip* 2ReA^- (26) 

Here ip is the semiclassical state and ///is the measure, a path 
integral over the fields in ip. // need not be known in detail, except 
that it enforces the reality condition in eq. (|26| ) ( via (— h)<5/i/<5E^ = 
2ReAf-/i). Also, /i must be real (//* = //) in order for norms to be 
real . If one inserts the relations eq. Q213| ) on the left in eq. ([24]) and 
carries out the complex conjugation (using fi* = //); the result is 
the right-hand side of eq. (24).D 



D Kinematic vs. physical operators 

I continue with the comparison of the BPR operators to field strengths 
in QED. In QED, operators representing observables should involve 
only true, rather than gauge degrees of freedom, therefore should 
commute weakly with the smeared Gauss constraint. Similarly, in 
quantum gravity, observables should commute with the constraints. 
The BPR operators do not do this. 

An operator will be said to be kinematic if it commutes only with 
the Gauss and diffeomorphism constraints; and physical if it com- 
mutes with all the constraints, including the scalar constraint. It is 
desirable for an operator to commute with as many of the Hamilto- 
nian constraints as possible, since the operator is then presumably 



14 



independent of one's choice of arbitrary coordinates. Physical oper- 
ators are best of all, but even a diffeomorphism and Gauss invariant 
operator can furnish valuable information about the physical mean- 
ing of a state. The volume, length, and area operators, for example, 
are only diffeomorphism and Gauss invariant. 

I shall show below that the BPR amplitudes for the unidirec- 
tional case can be modified so as to make them kinematical. I shall 
refer to the modified amplitudes as the unidirectional BPR ampli- 
tudes, since they continue to refer to excitations moving in a single 
direction. (That is, it is not necessary to combine a right-moving 
operator with a left-moving operator to get a kinematical opera- 
tor.) Further, in both the scattering and unidirectional cases, it is 
possible to construct additional combinations of the BPR operators 
which are actually physical, rather than kinematic only; but (except 
in the linearized limit) these physical combinations no longer refer 
to excitations moving in a single direction. 

From the remarks at the end of section I, in the planar symmetry 
case only three constraints Q survive gauge fixing: the generator of 
spatial diffeomorphisms along z, the Gauss constraint for internal 
rotations around Z, and the scalar constraint. Consider first C z , the 
diffeomorphism constraint. The BPR fields are local, scalar func- 
tions of (z, t), of density weight zero, and therefore do not commute 
with C z . Evidently the BPR fields have unacceptable dependence 
on the arbitrary coordinate label z. The simplest remedy is to get 
rid of z by integrating the BPR amplitude F, over z, following the 
prescription that observables in quantum gravity should be non- 
local; but the density weight of F must be unity rather than zero in 
order for integration to produce a diffeomorphism invariant quan- 
tity. One would like a density weight of unity for another reason: 
operators which are not density weight unity often turn out to be 
background- metric-dependent when regulated P?|. 



There are two ways to modify the BPR operators to make them 
density weight unity. (Either) replace 

A± -> E a ± ,A±; 
-AJ + 2ReA^ - E£,[-A* + 2ReAJ]; a = x,y only; 

(unidirectional) (27) 
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(or) replace 



-A^ + 2ReA B 



E|A*; 

Eb[ — A B + 2ReA^]; a,b = x,y only; 



(physical) 



(28) 



Classically, both these transformations are invertible. (Recall that 
the matrix is block diagonal, with one 2x2 subblock containing a 
= x,y, and A = X,Y or +,- only.) Both sets of new operators contain 
the same information as the old, therefore; and I can adopt either 
set as the new quantum BRS criterion (after integration over dz 
and inclusion of appropriate holonomies to secure gauge invariance; 
see below). The first set is easier to interpret, since each opera- 
tor contains only A+ or only A~ but not both, hence the name 
"unidirectional" : each operator refers to either left- or right-moving 
excitations, as do the original BPR operators. Each operator in the 
second set is a mix of left- and right-moving operators. The first set, 
though not gauge invariant, can be made kinematic by sandwiching 
between appropriate holonomies and intgrating over z. The second 
set needs only an integration over z to make it (not only kinemati- 
cal, but also) physical. I now verify the foregoing two statements in 
detail. 

The first set can be made Gauss invariant by sandwiching each 
operator between holonomies. Since the internal rotations around 
X and Y are fixed, and only rotations around Z remain, the gauge 
group is U(l) rather than SU(2). The irreducible representations 
are one-dimensional and labeled by a single integer or half-integer, 
the spin along Z. In eq. ( ^7|) some of the new unidirectional BPR 
fields have spin ± + ±' = 0, and are already Gauss invariant, while 
others have spin ± + ±' = ±2; only the latter fields are non-invariant 
and need to be sandwiched between holonomies, for example 




xM(z,z 1 )M(z 1 ,z r ) 
(unidirectional) , 



(29) 
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where M is a holonomy along the z axis, 

M(z 2 , Zl ) = exp[i r SzAZ(z')dz'. (30) 

J zl 



S±, S z are the usual Hermitean SU(2) generators. In eq. fl2"9|)the 
integration over z has been added to enforce spatial diffeomorphism 
invariance. If M(z2,zi) is pictured as a flux line running from zi to 
Z2- then eq. (|29| ) contains two parallel flux lines forming a hairpin 
contour, with the open end at the right-hand boundary z r : reading 
from left to right, each flux line starts at z r , runs to A+, then 
to zi, the closed end of the hairpin; then the flux line turns around 
and returns to z r . The flux line must be pictured as "open" at the 
z r end, because the S_S_ in the middle of the contour changes the 
eigenvalue of the S z in the holonomies. For the ingoing holonomy 
M(z r ,z), S z in eq. ( |3~0|) is evaluated at one eigenvalue, S z = m z ; 
while for the return holonomy M(z r , zi), S z is evaluated at a different 
eigenvalue S z = m z + 2. At the open end of the flux line, the m z 
values do not match, and there is no possibility of taking a trace. 
The operators of eq. ( p9[ ) are gauge invariant, despite the absence 



of a trace and the presence of open flux lines at z = z r . (This is the 
point in the argument where one uses the assumption of unidirec- 
tionality.) The wave is travelling toward z r but has not yet reached 
that point; therefore the boundary condition at z r is flat space. In- 
ternal Gauss rotations must respect this boundary condition, which 
means the triad at z r cannot be rotated. Technically, the Gauss 
constraint Cq must be smeared, / A(z)Cc(z)dz; and the smeared 
constraint commutes with the operators in eq. (^) at z = z r , be- 
cause Lim z ^ Zr A(z) = 0. 

It is straightforward to rewrite the holonomies as exponentials to 
the m z or m z + 2 power, and convince onself that the value of the 
constant m z does not matter. Also one can use any (2j+l)x(2j+l) 
dimensional representation of the matrices S_, S z , and the value of j 
does not matter: changing j merely changes the operators in eq. (|29|) 
by constants. 

One can also check that the criterion of eq. ( p8| ) is physical, after 
integration over dz. 

Gl = J dzE|A*. (31) 
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(Husain and Smolin [|2T| use the notation Kg for these operators, but 



it is perhaps better to reserve the letter K for extrinsic curvature.) 
For the most part the check is a straightforward computation of the 
commutator between the criterion and the scalar constraint; and the 
proof has been done already by Husain and Smolin. I will discuss 
only the one step which involves a slight subtlety, an integration 
by parts. (This step is trivial in the case studied by Husain and 
Smolin, since their z axis has the topology of a circle and there are 
no surface terms.) One starts from the integrated scalar constraint 
/dz / 5N(z , )C s (z') , where 

C s = E C C E* e cd e CD A«A»e-e MN /4 + e CD E c c F°% (32) 

e c d is the two dimensional constant Levi Civita symbol (c,d = x,y 
or X,Y depending on context), and F is the field strength 

F Z D C = (d z 6» - 6 DE A Z Z )A C E . (33) 

<5N is a small change in the densitized lapse. Since it must leave 
invariant the coordinate system at the boundaries, it vanishes at 
boundaries. In II-III I used a scalar constraint which equals eq. ( |3"2"| ) 
divided by a factor of E|, but here I use eq. ( |3"2D itself, for simplic- 
ity, because factors of E| do not matter: the operators of eq. (|28[) 
commute with E|. The canonical commutator is 

[A^,E^]=M^ a b 5(z-z'), (34) 
and one must use this to evaluate 

[|dzE^,|dz'5N(z')C s (z')]. (35) 



Cs, eq. (|3"2"|), is the sum of two terms, and only the commutator 
with the second term involves an integration by parts. The second 
term is proportional to the field strength F° , so that the commuta- 
tor has the form 

|dzdz'5N(z') •■■ [E|AB(z),E c c F°(z')] 

= f dzdz'5N(z') •■■ {E c c [%,FM + E-[A^,EyF°} 
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= I dzdz'5N(z') ••• fi{E£(zO[-d z ^z-zOA°( Z )+e DB AX] 

+E a c [5(z-z')]F°} 

= J dzdz'5N(z') •■■ h{E£(z')[-F^(z-z')] 

+E a c [5(z-z')]F°} + ST 
= + ST, (36) 

where ST denotes the surface term which arises on converting — d z '5(z— 
z') to +d z 5(z — z') and integrating by parts with respect to z. 

ST = =/dz'5N(zO---{E a c (zO[5(z-zOAD(z)]-^ 

= 0. (37) 

At a first glance one might think that surface terms vanish only when 
integrating by parts with respect to z', since the <5N depends on z', 
not z. Because of the S(z — z'), however, 5N(z') kills surface terms 
arising from integration by parts with respect to either varable, z or 
z'. □ 

The discussion just given proves that the first constraint J E|A^ 
in eq. (|28|) is physical, but says nothing about its complex conjugate 
partner /Efj[— A^ + 2ReA^]. In the previous section and in part B 
of the present section I remarked that the BPR constraints come in 
complex conjugate pairs, for example, 

A" = 0; 

-A+ + 2ReA+ = 0. (38) 

The second member of the pair is much harder to work with, be- 
cause of the non-polynomiality of the ReA^ factor. This , plus 
the complex conjugate pairing, suggests that an indirect approach 
might save labor: do the calculation for the first constraint, then 
take the Hermitean conjugate to obtain the corresponding result for 
the harder constraint. For example, commute the first constraint 
in eq. with Q, then take the Hermitean conjugate of the com- 
mutator to obtain the commutator of the second constraint with 
Q . This indirect approach works when Q is the Gauss or diffeo- 
morphism constraint, because these two constraints are self-adjoint. 
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However, the scalar constraint is not self-adjoint for my factor order- 
ing, or indeed, for any of the factor orderings most often encountered 
in practice when working with a complex connection. Therefore one 
cannot simply take a Hermitean conjugate to prove that the second 
constraint in eq. (p8|) is physical; one must directly work out the 
commutator of the second constraint with Cg. I now do this. 

Proof that the second constraint of eq. (ES|) is physical^ since the 
sum of the two constraints of eq. (p8|) is proportional to Eg Re A a , 
one need only prove that this expression (or more precisely its inte- 
gral over z) commutes with the scalar constraint. The Re A a factor 
is non-polynomial; but fortunately the calculation is not involved, 
because Eg Re A a turns out to be a total derivative. Therefore its 
integral is a surface term, and the surface term is physical, because 
Cs is smeared by a <5N which vanishes at the surfaces. The following 

formulas may be used to prove that E^ Re A® is a total derivative: 

ReA^ = e BJ ^ a ZJ ; 

u l J = [^i[ja] + %ai] - ^a[ij]]e lI e jJ ; 
% a] = e iM [«9jef -d. A ef}/2. (39) 

The first line is the reality condition for the transverse connections 
Ajj 5 . The Lorentz connection c<j a J is four-dimensional; e;M is the 
tetrad, and ef 1 is its inverse. The standard Lorentz gauge fixing 
condition, e tz = 0, may be used to simplify the sums over i in the 
definition of u; a J , when I = Z; one gets 

ujf = - ( 9 z g ij e zZ e jJ /2; 
E^ReA* = ee h B e BJ u z J 

= -<9 z g aj e bj /2, (40) 

which is a total derivative as required. □ 



It is fortunate that both constraints in eq. (|28 ) are physical; oth- 
erwise the scalar constraint would treat the two polarizations of 
gravitational waves differently, since the two constraints are needed 
in eq. (pg|) because there are two polarizations. It would not be sur- 
prising if the two polarizations are treated differently in the wave- 
functional, since that quantity is not directly observable. In fact the 
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two polarizations are treated differently, at least in the wavefunc- 
tional for the linearized limit J28] . However, it would be surprising if 
the two polarizations have different behaviour under commutation 
with the scalar constraint. That would lead one to question the 
usual choices of factor ordering, since this ordering causes the scalar 
constraint to be non-Hermitean, and a Hermitean constrain would 
treat the two polarizations symmetrically. . In section VI, I discuss 
indirect evidence that the factor ordering should be changed; but I 
have no direct evidence that the ordering is incorrect. 

(This paragraph probably can be skipped on a first reading.) 
In the plane wave case one has special (z,t) coordinates, call them 
(z',t'), which parameterize the hypersurfaces z' ± ct' = const. One 
can therefore envisage solving the problem of meaningless z coordi- 
nate in another way, by a transformation of coordinates to the more 
meaningful coordinate z'. One might still integrate over dz, but with 
a factor of dz' /dz inserted in the integrand. (This factor would give 
the desired density weight of one.) I have perhaps not given this 
option the attention it deserves. Note, however, it would require a 
great deal of attention. Information about the special surfaces z' = 
const, is encoded in the metric, so that the factor dz' /dz is not just 
a function, but a functional of the metric components ( and a highly 
non-polynomial function at that!) This factor becomes an operator 
in the quantum case, and it is a highly non-trivial question how one 
regulates and factor orders this operator. This is the problem one 
runs into if one shifts to z' at the quantum level. If one trys to shift 
at the classical level, then one must fix the C z gauge and replace 
Poisson brackets by Dirac brackets which respect the gauge-fixing 
condition. The classical Dirac brackets between the surviving ob- 
servables tend to be messy, and so far I have found no quantum 
representation. 

(The following three paragraphs contain material which at first 
glance may seem to be of only historical interest, but will be needed 
later in section V.) I interpret the BPR criteria in a semiclassical 
sense, as < F >= rather than F-0 = 0. Yet the Hamiltonian 
constraints are always imposed strongly, as C^p = 0, even though 
(in the linearized limit, at least) the Ci are sums of creation and 
destruction operators, like the BPR field strengths. Why this differ- 
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ence in treatment? This same question was posed and answered in 
a different context, Lorentz gauge QED, many years ago, and it is 
worthwhile to take a moment here to review that discussion 1201. In 



Lorentz gauge QED, the analog of the Q is (the usual Gauss con- 
straint, plus) the four- divergence dA = 0. The analog of the strong 
requirement Ci^ = would be dAip = 0; and the analog of the 
semiclassical requirement < Q >= would be dA + ip = 0, where 
the superscript + denotes positive frequency components. (Since a 
splitup into positive and negative frequencies is available in QED, 
there is no need to introduce a semiclassical average.) 

Both constraints, the strong and the positive frequency /semi classical, 
are used in the Lorentz gauge literature. Authors who employ the 
positive frequency constraint tend to treat the "unphysical" part of 
the Hilbert space with more respect. (Remember that the Lorentz 
gauge condition is designed to eliminate the effects of the unphysi- 
cal, longitudinal and timelike components.) Heitler [3(| is a typical 
proponent of this approach: he gives a very careful treatment of the 
unphysical sector, including a full discussion of the Gupta-Bleuler 
formalism. The payoff is that dot products over the full Hilbert 
space are well-defined, including dot products of longitudinal and 
timelike photons. Authors who employ the stronger constraint [[II 



pay a price: it is possible to find states which are annihilated by 
both the annihilation and the creation parts of dA, but these states 



are not normalizable in the unphysical sector |29|. This result is not 



particularly surprising: since the creation operators in dA create a 
state with one more timelike or longitudinal photon, ip must be a 
sum over an unbounded, infinite number of longitudinal and time- 
like occupation numbers. This infinte sum leads to the divergence 
in the norm. The authors who use the strong constraint are well 
aware of this difficulty, and they circumvent it by requiring that the 
dot product in Hilbert space be taken over physical excitations only. 

Returning to the gravitational case, one can now see why the 
strong criterion will work for the Q, but not for the BPR operators. 
For the moment, imagine the gravitational theory to be linearized, so 
that the analogy to QED is strongest. The creation operators for the 
BPR operators create physical quanta, not unphysical. If I impose 
the strong criterion, I get states which have an infinite norm in the 
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physical sector. There is no way of avoiding this by restricting the 
measure at a later step. If I now pass from the linearized to the full 
theory, there is no reason why the strong criterion should suddenly 
become applicable. I must use the semiclassical criterion, which is 
justified using the correspondence principle. 



IV Additional Solutions 



In the previous section I derived quantum BPR operators. In the 
present section I construct new solutions to the constraints. In the 
next section I apply the BPR, ADM energy, and Lz operators to the 
solutions constructed in paper II-III and this section. 

I start from the solutions considered in III. These are strings of 
transverse E A operators, ordered along the z axis, and separated by 
holonomies: 

-TT f Zn+1 

^ = Q.I / dz i e(z i+ i-z i )M(z i+ i,z i )E A i .(z i )S A i6(zi-zo)]M(zo,z n+ i). 

(41) 

The M are holonomies along z, 

M(z i+1 , Zi ) = exp[i f n+1 A z z (z')Szdz'], (42) 

J Zi 

and the Sm are the usual Hermitean SU(2) generators. These can 
be 2j+l dimensional; they need not be Pauli matrices. The func- 
tions in eq. (fl ) are Heaviside step functions which path-order the 



integrations, z < zi < • ■ ■ < z n+ i. For this section only, the bound- 
ary points zi and z r are relabeled zo and z n+ i . Although the metric 
is not flat at the boundaries, it can be taken as conformally flat at 
boundaries, with any radiation present confined to a wave packet 
near the origin ||. 

Since the full SU(2) invariance has been gauge fixed to 0(2), it 
is convenient to use basis fields introduced at eq. (^), fields which 
are irreducible representations of 0(2). These are one-dimensional, 
labeled by the eigenvalue of Sz, e. g. , 



= (E| ± iE Y )/V2; 
E A S A = E*S_ or E1S+. (43) 
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Because the irreducible representations are one-dimensional, there 
is no need to sum over both values of A ; = ± in eq. (fUD, in order to 
obtain a Gauss-invariant expression; nor is it necessary to take the 
trace in that equation. However, one must be sure to have an equal 
number of S + and S_ matrices in the chain, in order to form a closed 
loop of flux with no open ends violating Gauss invariance. That is, 
if one visualizes each holonomy M (z i+1 , z ; ) as a flux line along z from 



z; to z i+1 , then the factor in the square bracket, eq. ([|l]), may be 
visualized as a flux line from z to z n+1 . The line varies in thickness 
(varies in Sz eigenvalue) because of the S± operators encountered 
along the way, but the final Sz value at z n+ i must equal the initial Sz 
value at z (there must be an equal number of S + and S_ matrices 
in the chain). Then the final holonomy in eq. (ffTD, M(zo,z n+ i), can 
join the two ends at zq and z n+ i and turn the open flux line into a 



closed flux loop. As shown in III, the wavefunctional of eq. (II) can 
be made to satisfy all the constraints, by suitable choice of the ai 
and A;. 

In order to obtain a larger, kinematical space, as well as more 
physical solutions, one can simplify (!) eq. flUD by dropping all Q 
functions. This removes the path ordering, or ( in visual terms) this 
allows flux lines to double back on themselves. 



i=l 



Am = II f n+1 dz i M(z i+1 ,z i )E^(z i )S A iM(z ,z n+1 ). (44) 

zO 



Again, the expression is Gauss invariant even if A; = ± is not 
summed over; and no trace is needed. 

To check the spatial diffeomorphism and scalar constraints, one 
must first obtain these constraints from the Hamiltonian, written 
out in an 0(2) eigenbasis: 

H T = N'pl (EirW+A,- + E(±i)E b ± FS 

± 

+iN z ]TE^ b 
± 

-iN G [d z E z z-^±i)E^] + ST 

± 

= N'Cg + N Z C Z + N G C G + ST, (45) 
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where 

F z T b = % T iAj]Af. (46) 

( 2 ^E is the determinant of the 2x2 transverse subblock of the matrix 
E^. ST denotes surface terms (terms evaluated at the two end- 
points on the z axis, zo and z n+ i). The detailed form of these terms 
is worked out in II but will not be needed here. The primed lapse 
N' equals the usual lapse N multiplied by a factor of E|, and corre- 
spondingly the scalar constraint Cs is the usual constraint divided 
by E|. As shown in II, this renormalization leads to a much simpler 
constraint algebra, but again the details of this will not be relevant 
here. The system is quantized by replacing transverse A^ and E| by 
functional derivatives, these being the fields conjugate to the fields 
in ip. 

At - M/<% 

E z z - -M/5A Z Z . (47) 

The operator ordering (already adopted in eq. fl4"T|)) is functional 
derivatives to the right. The first term in eq. Q4"5p contains an inverse 
operator (E|) _1 ; this is well-defined provided E|M never vanishes, 
that is, provided the Sz in eq. ( |4"2| ) never has the eigenvalue zero. 

The physics must be invariant under small changes SN in the 
lapse and shift, so that the constraint should be written as 

= J dz^N'Cg + 5N z C z ]^kin; (48) 
= 5N'(z b ) = 5N z (z b ), (49) 

where z b is either boundary point, z or z n+1 . Eq. ([13) guaran- 
tees that the boundary conditions at z b are left unchanged by the 
transformation of eq. (pESp. 

Both the Cs and C z constraints in eq. (|48D contain terms propor- 
tional to F^, the field strength defined at eq. fl46|) . When a typical 
term of this type acts on ip, the result is (up to constants) 

J dz5N(z)E;F z - b [Vkin] 

J dz5N(z)E^F z - b [- ■ ■ J dziME^S.M ■ • •] 
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= J dz6Nj2^+[---J dziM(9 z -iA z z )5(z-Zi)S_M---] 

= ...[..•/ d Zi 5(z - Zi )(9 zi - iA z )MS_M- • •] 

= •■•[•■• J d Zi 5(z - Zi )(-iM[S z , S_]MA Z - iA z z MS_M) • ■ ■] 
0. (50) 

On the third line I have changed the <9 Z to <9 z j and integrated by parts 
with respect to z;. The surface terms at z ; = z b vanish because the 
5N(z)5(z — Zj) yields a factor of 5N(zb), which vanishes at bound- 
aries. Thus ip\a_ n is annihilated by all constraint terms containing 
field strengths F^. This is already enough to prove that the spatial 
diffeomorphism constraint C z annihilates ^kin, hence ?/>km is at least 
part of a kinematical basis, if not a physical basis. 

The state would be physical if it were also annihilated by the 
first term in Cs, which I call Ce because of the ( 2 ^E factor which 
it contains. (This term also happens to be proportional to a field 
strength, namely F z y .) The following is a sufficient condition for 
V^kin to be physical: Ce annihilates the state if it contains only two 
out of the four transverse fields: 

(either) E* and E* only ; 

(or) E+ and E y _ only. (51) 

If t/> contains E+ and E* only, for instance, their product E* E* 
is not contained in the triad determinant ^-E ; therefore the con- 
nection determinant e a bA+A^ in Ce will necessarily annihilate any 
wavefunctional containing only E* and E* . Note the wavefunctional 
must have an equal number of E* and El fields, because Gauss in- 
variance requires an equal number of S_ and S + operators in the 
chain. 

One can generate additional physical states, starting from those 



described by eq. fl44|) and eq. (pT|), by applying the operators G* 
or G£, where the G^ are integrals over z of the weight one objects 
constructed at eq. (|28|): 



n+l 

zO 



G^ = / dzE^. (52) 
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As shown in section III, the operators commute with H; they 
are physical. Hence application of m factors of G^ to a functional 
^[E*,E*] replaces m x superscripts in the chain by y superscripts, 
but leaves ip physical. The operators G^ change the total intrinsic 
spin of the wavefunctional; for a discussion of this point, see section 
VI. 

I have labeled the generators S^i in eq. (|l|) using two quantum 
numbers, j and m (more precisely, initial or final m, if the generator 
SaIs one which changes m). Once the SU(2) symmetry is broken 
to 0(2), however, the j quantum number loses signifigance. I could 
replace the Sai by any other matrix with the same m (and Am) but 
different j, and ip would change only by a constant factor. 

Even though the j has no physical signifigance, it is mathemati- 
cally convenient to use Sa having definite j. One can then employ 
the familiar commutation relations of the Sa in calculations. Also, 
the planar state ip is presumably a limit of some three-dimensional 
state for which the label j has meaning. The fact that states of 
different j are equivalent in the planar limit presumably means that 
the correspondence between three-dimensional and planar states is 
many to one. 



V Application to Solutions 

In this section I study the solutions constructed so far by applying 
several operators to them: the modified BPR operators (from sec- 
tion III), the ADM energy operator (from paper III and appendix 
C), the area operator E| for areas in the xy plane (from paper III), 
and the operator Lz giving the total spin angular momentum around 
the z axis (from paper IV and appendix B). All the solutions (those 
constructed in papers II-III as well as the new solutions constructed 
in section IV) have the form of strings of transverse E^z^S^ oper- 
ators separated by holonomies M(zi + i,Zi) . The solutions in papers 
II-III have additional step functions 9(z 1+ \ — z;) which path order 
the integrations over the z ; , but the 9 factors will play a minor role 
in the considerations of the present section. 

Consider first the ADM energy operator. Often this is identi- 
fied with the surface term in the Hamiltonian, but, as discussed in 



27 



appendix C, the volume term can also contribute. In the present 
case the volume term typically does contribute, but its only effect 
is to double the size of the surface term, and I will ignore volume 
contributions. The surface term is, from eq. ( p4j) , 

Hst = — ^MNEM^b Izi 

= ih[E h _8/5E h _-E\8/8E\}^ (53) 

When this operator acts upon a factor of E^_(zi)dzi in the wavefunc- 
tional, it gives ^ihE^dz; times a factor of 5{z\ — z r ) or <5(zj — z\). Ob- 
viously none of the solutions is an eigenfunction of the ADM energy, 
since the 5 function deletes one integration dz ; . One could perhaps 
construct an eigenfunction by summing over an infinite number of 
solutions, each containing one more dzi integration. Each additional 
integration should be multiplied by an additional factor of i, to can- 
cel the i in eq. and make the eigenvalue real. Investigation 
of such sums is beyond the scope of the present paper. Without 
a measure one does not know whether such a sum converges to a 
normalizable result. 

If the Gauss constraint, 

C G = -i[«9 z E z z -e M NE^AN], (54) 

vanishes at the boundaries, it can be used to simplify the ADM 
surface term to 

fpz |z r 

st — — -t^Ziz U 

= h(5/5A z z ), z . (55) 

Unwanted factors of i are now more of a problem. A z occurs only in 
holonomies, where it is always multiplied by (a real matrix Sz times) 
a factor of i, and the E\ will bring down this factor of i. Except for 
the solutions constructed in paper II, there is always at least one 
boundary where the Gauss constraint holds, so that factors of i will 
be a generic problem. Of course one could eliminate the problem by 
discarding the holonomy structure, but this is a solution almost as 
unattractive as the problem. 

The operator E%, z occurs in the ADM energy, while E\ itself is 
the area operator for areas in the xy plane. (E| = ee| = ^e.) 
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Therefore the area operator also has pure imaginary eigenvalues, a 
situation already noted in Appendix D of paper II; see also DePietri 
and Rovelli [0. 

Even if one for the moment ignores the factors of i, there is an- 
other problem with the area operator: at any boundary where the 
Gauss constraint is satisfied, the area operator will always give zero. 
If the Gauss constraint is satisfied, say at the left boundary zi, then 
there is no net flux exiting at z\. One can regroup the holonomies 
until there are no M(zi, z{), only M(z r , Zi); or until every M(z ; , zi) 
is paired with an M(zi, zj) to give M(z;, zj M(zi, zj) = M(z;, zj). 
Either way, there is no holonomy depending on (zi), and the area 
operator E|(zi) gives zero. 

Next consider the action of the BPR operators. The solutions 
given in II (and some of those considered in section IV) contain 
either E+ operators or Et operators, but not both. A wavefunctional 
which contains only E 3 ^ operators (for example) will be annihilated 
by A~ = M/5E+, even before any semiclassical average is taken: 

A-^Et) = 0. (56) 

In classical theory, A^ = is a signal that the solution is pure left- 
moving; and from this one might expect that ^[Ei] is unidirectional. 



Condition eq. fl56|) is too strong, however. As discussed at eq. (j2lD 
of section III, One expects at most a vanishing semiclassical average 
(A~) = 0. In fact from the remarks on Lorentz gauge QED in 
the concluding paragraphs of section III, eq. fl56|) implies that if) is 
probably not a normalizable state! 

The solutions of paper III and most of those from section IV 
contain both E^_ and Ei, hence are not annihilated by any BPR op- 
erator. One cannot conclude that these solutions are infinite norm, 
therefore. However, they do suffer from the problems described pre- 
viously in this section, those associated with the ADM energy and 
area operators. 

VI Conclusions 

Papers II-III proposed new solutions to the constraints, and section 
IV of the present paper proposes still more solutions. However, the 
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investigations of section V have demonstrated that these solutions 
are less than satisfying in several respects. 

This outcome is perhaps not surprising. In earlier work , sim- 
plicity was the primary criterion for choosing the method of quan- 
tization (polarization and factor ordering) , as well as the primary 
criterion for constructing solutions. Simplicity is the criterion one 
uses when information is scarce, however. In the present paper sev- 
eral operators of physical signifigance were available, and the theory 
and its solutions can be held to a standard more demanding than 
simplicity, the standard of a reasonable physical interpretation. As 
a result, solutions are called into doubt; but this happens for a rea- 
son which is fundamentally positive: more is known about how to 
interpret and understand the solutions. 

The difficulties with imaginary eigenvalues encountered in section 
V seem to require fundamental revisions in the theory, since the dif- 
ficulties are closely linked to the complex nature of the connection. 
As discussed in section V, the operator E|, present in both the area 
operator and the ADM energy operator, has complex eigenvalues 
because of the i in the holonomies, exp[iA^ • • •]. Getting rid of the i 
entails dropping or modifying the holonomic structure, not a pleas- 
ant prospect. Recently Thiemann [2(J has proposed an alternative 
formalism based upon a real connection. Thiemann's alternative is 
motivated primarily by issues of regularization, but has the desirable 
side effect of producing real eigenvalues for the area operator. 

It is a little harder to see how switching to a real connection 
will cure the problem of zero eigenvalues of the area operator at 
boundaries, also uncovered in section V. The zero eigenvalues occur 
only at boundaries where Gauss invariance is satisfied . Since E| 
is a Gauss invariant, at first sight any connection between area and 
Gauss invariance seems strange. The connection is indirect, via the 
structure of the complex connection A^ = i Im A^ + Re A^ . Im A^ 
is the part which does not commute with E|; therefore its presence 
in the wavefunctional gives rise to non zero area. Re A^ is the part 
which transforms like a connection under Gauss rotations; therefore 
it is needed in the wavefunctional for gauge invariance. The notions 
of area and gauge invariance are linked only because Im A and 
Re A are linked together to form a single (complex) connection. At 
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a boundary where Gauss invariance is satisfied, if there is no net 
flux exiting, then there is no dependence on the connection, hence 
no area. In the Thiemann scheme one still joins Im A and Re A 
together to form a single (real) connection; in fact the Thiemann 
connection is just the Ashtekar connection without the factor of i. 
The real and imaginary parts of the connection are separated when 
constructing the Thiemann constraints. Must they be separated 
when constructing the wavefunctional as well? If the answer is yes, 
the wavefunctional would not be purely a product of holonomies. 

Before doing anything as drastic as dropping the holonomic struc- 
ture, it is a good idea to investigate what happens to the zero area 
argument when it is extrapolated from the planar case to the full, 
three-space-dimensional case. The argument that Gauss invariance 
leads to zero area depends on properties of the wavefunctional at 
boundaries, and the behaviour at boundaries changes markedly with 
spatial dimension. 

In the full three-dimensional case, the smearing function for the 
Gauss constraint must vanish at spatial infinity, so that there is no 
need for the Gauss constraint to annihilate the wavefunctional there. 
As a result, net flux may pass through the boundary at infinity, and 
there is no difficulty obtaining finite area at the boundary, even 
when the wavefunctional is purely a product of holonomies. This 
suggests that the zero area may be a problem which occurs only in 
the planar limit. 

In fact the problem may not exist even in the planar limit, if one 
takes this limit correctly. Imagine the generic three dimensional flux 
configuration which is well approximated by planar symmetry: near 
the origin, the flux lines corresponding to holonomies containing 
dx a are finite in cross section and well-collimated along the z axis. 
The planar wavefunctionals constructed in papers 1 1- 1 1 1 contain fac- 
tors of S±, presumably relics of Clebsch-Gordan coefficients coupling 
A^ dx a holonomies to holonomies containing A^ dx a and Aj dx a . 
The latter are represented by flux lines lying in planes z = constant. 
Although A^ dx a flux lines are well collimated near the origin, they 
must diverge far out along z into the past or future. A sketch of 
the A^ dx a flux lines resembles a drawing depicting radial geodesies 
near a wormhole: the flux lines come in from radial infinity, pass 
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through a narrow "throat" oriented along z, then diverge once more 
to radial infinity. (The wavefronts perpendicular to these rays are 
constructed from dx a and A J dx a holonomies.) Alternatively, 
the A^ dx a flux lines at infinity may not exit through the surface at 
infinity but may loop back and close on themselves, resembling the 
flux lines of a solenoid in magnetostatics. Either outcome is allowed 
by the boundary conditions on the Gauss smearing function at in- 
finity, and for either behavior at infinity, the behavior at the throat 
is the same. If one takes a cross section through two points z\ and 
z r > z\ at the throat, one finds net A^ dx a flux through both bound- 
ary points z\ and z r . If this picture of the three dimensional flux is 
correct, then in the planar limit one should not impose Gauss invari- 
ance at the boundaries. Planar solutions would resemble the "open 
flux" solutions studied in paper II. If there is A^ dx a flux through 
the boundaries, the zero area problem at boundaries disappears. 

Even though the A^ dx a flux lines now extend throughout the 
entire range z\ < z < z r , one can still construct localized wave pack- 
ets. In paper II I reviewed the geometrodynamical treatment of the 
planar problem, and introduced the Szekeres scalar fields B, W, and 
A. (It is a little easier to work out the boundary conditions for B, W, 
and A, rather than work directly with the Ashtekar fields; as shown 
in II, the boundary conditions on B, W, and A then imply corre- 
sponding boundary conditions on the Ashtekar variables.) When 
only the field A is present, the forces on a cloud of test particles are 
isotropic in the transverse (xy) direction; the elliptical distortions 
characteristic of gravitational waves appear only when the B and 
W fields are non-zero. (In a more covariant language, the compo- 
nents of the Weyl tensor which give rise to transverse deviations of 
geodesies are present only when B and W are non-zero.) One can 
impose wave packet boundary conditions on B and W (equivalent ly, 
on transverse components of the Weyl tensor), requiring these quan- 
tities to vanish at the boundaries z b ; but this requirement tells us 
nothing about the behavior of A at the boundaries. The variable A 
determines geometrical quantities such as areas: the Ashtekar area 
operator E| is just exp(A). Perhaps one has a "wave packet", but 
not a "geometry packet". One may require localized, wave packet 
behavior for B and W, but not for the more geometrical quantity A. 
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To summarize, there are two possible solutions to the zero area 
problem. The first splits the connection, abandoning the strict holo- 
nomic form for the wavefunctional. The second allows the Gauss 
constraint to be non-zero at boundaries and in effect assumes that 
the "open flux" boundary conditions used in paper II are generic. 
Further information and thought is meeded before one can decide 
between these two alternatives. 

Apart from the zero area difficulty, there is another reason why 
the transition from three to one space dimension needs more atten- 
tion. Ultimately one would like to use the planar case as a guide to 
the behavior of radiation in the full, three dimensional case. In the 
full case, one expects the connections to occur in holonomies, so as 
to preserve gauge invariance. In the planar case, the gauge fixing 
allows the transverse connections to occur outside of holonomies. 
Two key radiative properties of the planar solutions, their direc- 
tionality and spin, are associated with the transverse sector, which 
least resembles the three-dimensional case. It will probably be nec- 
essary to recast the transverse sector in a more holonomic language, 
in order to understand more clearly what happens on passing to the 
full theory. 

For the moment let us overlook any possible difficulties with zero 
area and suppose that one shifts to a real connection, in order to 
eliminate the problem with imaginary eigenvalues. One can ask 
whether the solutions constructed in papers II-III and section IV 
are likely to survive the shift to a real connection formalism. The 
present solutions may not survive, if the factor ordering is changed; 
and it is easy to imagine a reason why one might want to change the 
factor ordering. The scalar constraint usually must be taken to be 
non-Hermitean, in a complex connection formalism; whereas with 
a real connection one may wish to factor order so as to make the 
constraint Hermitean. 

It may be helpful to comment briefly on why the scalar constraint 
is difficult to make Hermitean in a complex connection formalism. 
The usual recipe for constructing a self-adjoint operator is to factor 
order it, and then, if the operator is not self-adjoint, form the aver- 
age (C + Ct)/2. Cf is C, with the order of all operators reversed and 
the connections A replaced by At; in turn the At are replaced by -A 
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+ 2 Re A. This last step introduces the unwanted non-polynomial 
expressions Re A into the Cf term. In the case of the Gauss and spa- 
tial diffeomorphism constraints, identities may be used to eliminate 
the Re A contributions, and Cf is well-behaved, in fact identical to 
C. In the case of the scalar constraint, the unwanted Re A terms do 
not go away. Within a real connection framework, the Af is just A, 
and the traditional (C + Cf)/2 recipe is easier to implement. 

Although the present solutions may not survive as exact solu- 
tions, they may constitute approximate solutions to the new, modi- 
fied constraints, solutions valid in the limit h— »0. This would happen 
because the new and old scalar constraints presumably will differ 
only by a reordering of factors, hence will differ by terms of order 
h. Also, the Gf, operators, defined at eq. (|3lD , wer shown to be con- 



stants of the motion by Husain and Smolin [^T|. These operators 



are likely to remain constants of the motion, in any transition to a 
new operator ordering, because of the close connection between the 
G£ and total spin |DJ. 

Even though the complex connection formalism may not be ap- 
propriate for the dynamics, this formalism is the natural one to use 
when constructing a criterion for the presence of radiation. Note 
that the BPR operators were derived in section II using only sym- 
metry considerations; no assumptions were made about the factor 
ordering or dynamics. Even if one dropped the Ashtekar connection 
and used the Thiemann connection, one would have to reintroduce 
the Ashtekar connection in order to express the results of section II 
succinctly! Anyone familiar with the classical results on radiative 
criteria will not be surprised at this: much of that work is most 
conveniently expressed using the language of complex connections. 
(See for example the work on the Weyl tensor quoted in section I 
and Appendix D.) 

My original intention in constructing the BPR operators was to 
elucidate the spin and directional character (that is, left vs. right- 
moving character) of radiative quantum states. Since these opera- 
tors turned out to be semiclassical in character, and no measure is 
yet available, it is perhaps early days to make a final judjment on the 
efficacy of the BPR operators. However, it is probably true that any 
unidirectional operator (BPR or other) will be at most kinematical; 
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that is, it will commute with the Gauss and diffeomorphism con- 
straints, but not the scalar constraint. Conversely, any operator 
which is physical (commutes with all constraints) will not be unidi- 
rectional; i. e. it will mix BPR operators of opposite directionality. 
The concept of unidirectionality is well defined in classical theory; 
see section II and appendix B. In quantum theory, however, every 
unidirectional wave must traverse a vacuum filled with zero point 
fluctuations which are moving in both directions. Presumably this 
is the intuitive reason why the modified BPR operators constructed 
in section III are either unidirectional or physical, but not both. 
Consistent with this interpretation of unidirectionality as primarily 
a classical concept, it is possible to prove that the unidirectional 
BPR operators are physical in classical theory; but factor order- 
ing problems prevent the proof from going through in the quantum 
case. Presumably any criterion for unidirectionality will have to be 
at most a semiclassical one. 

Classically, exact plane wave solutions are known in which the 
area operator E| = ^ 2 ^e evolves to zero [[RJ. In fact this collapse 



behavior appears to be generic; solutions which do not collapse are 
rare and are unstable under small perturbations |53]. The zero area 
cannot be removed by a change of coordinates, since typically there 
is an accompanying singularity in a scalar polynomial quadratic in 
components of the Weyl curvature tensor. One can ask whether 
a quantum- mechanical effect might prevent this collapse. It is not 
possible to answer this question definitively within the present con- 
text, because the area operator has imaginary eigenvalues. Nev- 
ertheless, one can see the outlines of a possible quantum solution 
which would avoid a collapse. The quantum area operator E| acts 
on holonomies exp(i/ Sz dz); so long as Sz is not allowed to 
assume the value zero, E| cannot have the eigenvalue zero. In the 
solutions constructed in papers II-III and section IV, the Sz value 
in each holonomy does not evolve dynamically, therefore remains 
non-zero if chosen to be non-zero initially. It remains to be seen 
whether this happy state of affairs will persist to a new formalism 
with real eigenvalues for the area operator. 
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A Details of the BPR calculation 



This appendix solves eqs. (|TlD and fll2|) for the connection and 
tetrads obeying BPR symmetry, the invariance group for unidi- 
rectional plane gravitational waves. When setting up a complex 
connection formalism, it is necessary to choose three phases: when 
defining the Lagrangian at the four- dimensional level, one must 
choose the duality phase 5 and the phase of exxYZ (see for exam- 
ple eq. (p~3|) ) ; and an additional phase comes in when rewriting the 
four-dimensional formalism in 3+1 canonical form.. These phases 
are explained in appendix A of II, and I use the same phase choices 
here as in that paper. I Begin at the four-dimensional level by solv- 
ing eq. (|TT|) for the Lorentz transformation L. 

Lpi = -d P eeave? . (Al) 

I have dropped a £ x d\ term which is zero because eq. (^), is a 
linear combination of d x , d y , and <9 V , all of which annihilate ef . In 
Rosen gauge, from eq. (|), 

dpC = SX a - ^g ua . (A2) 

Therefore 

L n = -e^ej 1 + e^. (A3) 

L is antisymmetric, as it should be. Eq. (|TTD determines L, eq. (|A3| ) , 
but otherwise imposes no new constraints on the tetrads beyond 
those already imposed by d x , d y , and <9 V = 0. 

Next consider eq. (0). In principle, one should be able to de- 
termine all the constraints on the A^ by solving these equations 
directly, but they are awkward, and it is easier to adopt an indirect 
approach. Given the tetrads, compute the Lorentz connection u^; 
then compute ^A , which is just the self-dual version of lo. In this 
way one finds that many components of ^A are identically zero. 
When this information is inserted into eq. flT2|), that equation re- 
duces to the trivial statement = 0, for most values of the indices; 
for a small number of index values the equation is non-trivial and 
can be solved with moderate effort. 

The equation relating u to the tetrads is 

^ja = eiiejja£ J 
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= [-gaj,i + gai, j + e j K«9 a e| < ]/2. (A4) 

From this equation, at least one of i, j, or a must be u, since deriva- 
tives with respect to x,y,v are zero. Further, the tetrad matrix in 
Rosen gauge is 2x2 block diagonal, with the zt to ZT block contain- 
ing constants only. This implies that at most one of i, j, or a must 
be u. After stripping off the (block diagonal) tetrads, one finds that 
the only non-zero to are 

^Va^, (A5) 

where a = x,y only and A = X,Y only. Now compute ^ 4 ^A from 
u, using 

2 (4) Ai J = ^ + i5( e IJ MN /2e TXYZ R MN , (A6) 

where the duality eigenvalue tfexxYZ equals +1, given my conven- 
tions 5 = exxYZ == 1- The only non-zero 

(4) A 

components are 

( 4 )Af7 4) Af, (4) Ar + - (A7) 

At first glance one might think this list is too short; there should be 
non-zero ^A^ ± as well, because the e IJ MN term in eq. (|A6|) will map 
the VA indices on ujJ a into UB (B = X,Y or ±). Duality maps VX 
into VY, not UY, however, because of the identities 

ctxyz = cuvxy = c V vy> (A8) 

etc. , where the last, mixed index tensor with two V's is the one 
which occurs in duality relations. This explains why there are no 

(4) A U± 

; to see what happened to ^A a , note 

2^ = -a V± + W(^N/26 T XYz)-a MN 

= c a v± + i<He v v y e TX YZ V a V± 

= ^Hl±S). (A9) 

For my phase choice 5 = +1, ^A}[~ vanishes. 

From eq. (A~7), there is no need to consider a = v in eq. flT2]). 



I consider first a = a = x, y. The first term in eq. (|T2|) is a linear 
combination of <9 X , d y , and d v , which vanishes for any a. The second 
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term vanishes from eq. ( |A2| ) and the absence of any A = v component 
of ^A. Then eq. (|I~2| ) collapses to 

= + + ^, (4) A a ' J + C J r (4) A a J ' - 0. (AlO) 

From eqs. ( |A3| ) and (0), the only non-zero elements of L are 

£ ux = i£ UY = (e c x + ie Y )/2 = e c + /V2; (All) 

or 

C v+ = e c + ;£u- = 0. (A12) 

Inserting this and eq. ( |A7| ) into eq. ( |A10| ), one finds that all the 
a = a = x, y equations are trivially = 0. 

Finally, consider a = u. The first term in eq. ( |T0"D vanishes as 
before. The only IJ index pair which does not give = is IJ = 
VA, A = X,Y only, which gives 

= o + d u e «A VA + C y v (4) AJ| A + C% ^A vy - d u C VA - (A13) 



Use eq. ( |A2| ) to simplify the first term; use eq. ( |A12|) to simplify 



the remaining terms and to show that the A = - equation is trivial. 
Then 

= + g ca A^ + - 2e;A-+ + 8 u e c + 

= g ca A a / +_2e c B AB+ + <9 u e + . (A14) 

I have used the duality relation A vu = — iA XY , and A~ + = iA XY . 
On the second line, recall that a "plus' index always pairs with a 
"minus' index to form the two-dimensional dot product: egA B+ = 
e^_A~ + + 0. The A u field in eq. ( A14 ) is a linear combination of 



A z and A t fields, and the A t fields are non-dynamical Lagrange 
multipliers for the Gauss constraints. I therefore eliminate the A u 
field in order to obtain a constraint on the dynamical field A v+ . 
From duality and eqs. (IA1) and (|A5|) for uj, 



2ei 4 B } A u B + = e aB K B+ + i(5/2 e TXYz)4> u MN ] 
= [e i+ o; aju + 0] 

= e^tejK^e^/2. (A15) 
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I solve eq. ( |A14j ) for A VA and insert eq. (|A15|) . 

«A?+ = ^[e }K d u ef]/2-g ca d u el 

= eJ + [2e jK <9 u e a - d u g ]a }/2 - <9 u e a+ + <9 u g ca e+ 

= <9 ugca e c + /2. (A16) 

The right-hand side of eq. (|A16|) is proportional to the part of A^ 4 " 
which contains no iS factor. 



2-Re'^A^ (- a VX + kO 
= e iV (e Xj + ie 
= e +j <9 ugaj /2. (A17) 



Therefore 

= _(4) A v+ + 2"Re" {4) A^+ 

= (4) Ay+, for 5 = -1. (A18) 

The second line means that the first line is the four dimensional 
connection computed with the opposite choice for the duality eigen- 
value, 5 = — 1 rather than 5 = +1. 

This result is very easy to transform from Rosen to a general 
gauge in the z,t sector, since the "minus" and "a" indices in the x,y 
sector remain invariant under such a transformation. In order to 
maintain the gauge conditions eqs. @ and (]T0|) on the tetrads, it 
is mecessary to combine any four- dimensional diffeomorphism (z,t) 
— > (z'f) with a Lorentz transformation, as at eq. (|ll]). A short 
calculation shows that a very simple Lorentz transformation L\ = 
dt'/dz will maintain all the gauge conditions of eqs. (|D and (|T0D. 
For transforming (^A one needs the corresponding £; the only non- 
zero matrix elements will be C T Z and C\; or equivalently C^j, £V\, 
and Thus the coordinate transformation to the general gauge 
amounts to a Lorentz transformation which multiplies eq. ( |A18| ) by 
an overall factor. 

= £ V v £ + - (4) Ar + , for 5 = -1, (A19) 

where every C and every ^ 4 ^A is to be calculated using the 5 — - 
1 convention. From eq. ( A19 ), the quantity in eq. ( |A18| ) vanishes 
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in every gauge. Similarly, from eq. (|A7|), the following quantities 
vanish in every gauge: 

= -(^+ + 2 "Re" WA"+ 

= (4) Ay- 

= (4) A"-. (A20) 

The eqs. ( A20|) and ( |A18| ) may be rewritten as 

= ( 4 >AJ+ - 2"Re' (4) AJ+ 

= < 4 >Af+ - 2"Re' ^Af+ 

= ( 4 )a t ~ 

= (4) Af-, (A21) 



where every connection in eq. (|A21|) is evaluated using the 5 



+1 convention. For the opposite duality convention, 5 = — 1, ex- 
change (+*-►—) everywhere in eq. ( |A21| ). For 5 = +1 but left- 
moving rather than right moving waves, again exchange (+<->—) in 
eq. (|A23). 

So far the calculation has been carried out entirely at the four- 
dimensional level. The four-dimensional connection ^A is related 
to the usual 3+1 connection A by the following equation from Ap- 
pendix A of II. 

A* = -eSUr (A22) 
Then the four-dimensional eq. (|A21|) implies the 3+1 dimensional 
equations 

= A"; 

= -A+ + 2"Re"A+. (A23) 

Again, exchange (+•*->•—) for the opposite duality convention or left- 
moving waves. 

B Kinematics of the fields: spin 

From paper IV , the integral 

L z = ij dz[E[(A* - ReA{) - (x <-> y)] (Bl) 
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gives the total spin angular momentum of the wave, and is a con- 



stant of the motion fI0 |. The integral is over the entire wave packet, 
that is from z\ to z r . As in II, the fields and Weyl tensor compo- 
nents which produce transverse displacements of test particles are 
assumed to vanish at the boundaries, with support only in the region 
z\ < z < z r . It is at first sight surprising that any conserved quantity 
associated with the Lorentz group should be given by a volume inte- 
gral (integral over z) rather than by a surface term (term evaluated 
at the endpoints z\ and z r ). However, in the one-dimensional planar 
case, the extensive gauge-fixing in the x,y plane removes all gauge 
freedom, except for rigid rotations around z, and the x,y sector of 
the theory resembles special relativity rather than general relativity. 

This appendix rewrites the integrand of Lz in terms of the unidi- 
rectional fields defined at eq. fl27|), in order to understand the spin 
content of the latter. Introduce triads e^ and inverse triads e^, and 
write the integrand of Lz as 

E[ImA x -(x^y) = [(e^e Kx - (x<-y)]E£ImA;[ 

= [(ef jeK]x + e[jeK) x - (x~y)]E£ ImA^B2) 

On the last line the term antisymmetric in J,K is proportional to 
E^ Im A^ 6jk- This expression is part of the Gauss constraint 
<9 Z E Z + E^A^ejK = 0, which implies E^ImA^e JK = 0. Hence the 
term antisymmetric in J,K can be dropped. The term symmetric 
in J,K can be expanded in 0(2) eigenstates, keeping in mind that 
every + index must be contracted with a - index. The J ^ K terms 
are proportional to 

x "I" e — e +x = e B e Bx 

= *l 

= 0. (B3) 

Hence these terms can be dropped also. The surviving terms are 
products of tensors with J = K, therefore helicity ±2 in the local 
Lorentz frame, a reassuring result. 

L z = - J dz[e^e +x Et Im A~ + e^e_ x E* Im A+ - (x«->y)] 

= i J dzKe +x E« [A" + (A" - 2 "Re" A~) 
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+e^e_ x E a + [A+ + (A+ - 2 "Re" A+)} - (x«-y). (B4) 

On the last line I have written Lz in terms of the unidirectional BPR 
fields introduced at eq. ( P7| ) This expression for gravitational spin 
angular momentum, possesses the same coordinate times momen- 
tum structure as the corresponding expression for electromagnetic 
spin angular momentum: 

L em = (l/4vr) J d 3 x[E x A 

= - J d 3 x[li x A] (B5) 

That is , one can interpret the unidirectional quantities E A and 
E (A - 2 "Re" A) in eq. (|B"4j ) as momenta associated with waves of 
definite helicity. This parallel with QED does not extend too far: 
these "momenta' have nothing like free-field commutation relations 
with each other, or with the triad "coordinates" . 

It is now clear why one wants the two combinations E a A+ and 
(-A~ + 2 "Re" A~) El to vanish: these two constraints remove 
left-moving helicity ± 2 contributions from Lz. Why must the re- 
maining, helicity zero combinations vanish? The helicity zero com- 
binations are E1A+ and (-A~ + 2 "Re" A~) E a . They are complex 
conjugates of each other, so that by adding and subtracting them 
from each other one gets pure imaginary and pure real constraints 

= E|(A* - ReA*) - ie AB ReA£E|; (B6) 
= -i eAB (A^-ReA^)E| + E|ReAf. (B7) 

Now consider the classical equation of motion 

= -iE z Zt - 5H/5A Z Z 

= -iE z Zt -E B A*. (B8) 

On the second line I use the Hamiltonian of eq. fl45|) . I also use 
the unidirectionality assumption (for the first time in this section; 
Lz is the spin operator also for the scattering case) and evaluate 
the Hamiltonian in Rosen (or at least conformally flat) gauge. The 
metric in a general gauge has the form 

ds 2 = [(-(N') 2 + (N z ) 2 )dt 2 + 2N z dzdt + dz 2 ]g zz + x,y sector, (B9) 
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so that to obtain conformal gauge, one must take N z = and N' = 1 
in the Hamiltonian, where N z is the shift and N' is the renormalized 
lapse defined following eq. From the real part of eq. ( |B8| ), the 
E Re A term in eq. ( |B7| ) vanishes. The rest of this equation is just 
the imaginary part of the Gauss constraint, 9 Z E| + cabA^ E| = 0, 
and vanishes also. This leaves eq. flBl]|). The e Re A E term may be 
simplified using the real part of the Gauss constraint; the E (A - Re 
A) term may be simplified using the imaginary part of the equation 
of motion, eq. (|B8| ). The result is simply 

= -i{8 t + «9 Z )E Z Z , (BIO) 

in any conformally flat gauge. This equation is discussed further at 
eq. fllTD of section II. 



C The ADM energy 

It is a worthwhile exercise to express the ADM energy in terms 
of BPR operators. In the usual three-space dimensional case, the 
Hamiltonian expressed in terms of the original ADM variables is the 
sum of a volume integral plus a surface term H st El, KB] , 



H = J d 3 x[NC sc + N'Q] + H st . (CI) 

In the classical theory, the ADM energy is just the surface term, 
since the constraints in the volume term must vanish everywhere 
when the solution obeys the classical equations of motion. Often 
one says that in the quantum case the ADM energy is just the 
surface term also, but this is not quite right, as we shall see in a 
minute. 

In the planar, one-space dimensional case, the expression for the 
Hamiltonian in terms of Ashtekar variables looks superficially much 
the same as eq. ( |CT| ) M, 

J dz[N'C sc + N Z C Z + N G C G ] + H^, (C2) 

except for the additional Gauss constraint, and the prime on N'. 
(The prime means I have renormalized the usual Ashtekar lapse 
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by absorbing a factor of E| into the lapse, as explained in II.) In 
both one and three spatial dimensions, one might be tempted to 
drop the volume terms, in the quantum mechanical case, because 
the constraints Q are required to annihilate the wavefunctional. 
However, the statement that the scalar constraint (say) annihilates 
the wavefunctional means, not C sc ip = 0, or even / dzC sc ip = 0, but 
rather 

J dz5N'C sc ^ = 0, (C3) 

where SN' is a small change in the lapse. The arbitrary change SN' 
must preserve the boundary condition at spatial infinity, N'— >1 . 
Hence <5N'— >0 there. On the other hand, when C sc occurs in the 
Hamiltonian of eq. (|C2|), it is multiplied by N' rather than SN'. 
There is no need for N' to vanish at the boundaries (in fact it be- 
comes unity there). Now suppose the evaluation of the action of C sc 
on if) reqires an integration by parts with respect to z. In eq. (|C3|) , 



when the constraint acts upon ip, surface terms at z = Zb will vanish 
because of the boundary condition on SN'. When H acts upon ip, 
however, the C sc in H is smeared by N' rather than SN'; the former 
does not vanish at boundaries. Consequently the volume term can 
contribute to the ADM energy in the quantum case. In the planar 
case both the Gauss and scalar constraints in the volume term can 
contribute to the ADM energy; neither N' nor Nq is required to 
vanish at the boundaries. In the usual 3+1 dimensional case with 
flat space boundary conditions at infinity, only the scalar constraint 
in the volume term can contribute; the remaining constraints are 
smeared by N; which are required to vanish at spatial infinity. 

For the plane wave case, the surface terms in the Hamiltonian 
were computed in section 4 of II. 

H st = — £mnEmA" |^ 
= I[F> A+ - E% Ag] 

= i^E h _S/SE b _-E h + S/SE h + ]ll (C4) 

To simplify the boundary term quoted in II, I have invoked the 
boundary conditions N z — >0, N'-^l on the shift and renormalized 
lapse. Evidently the ADM energy contains the BPR operators which 
are sensitive to the long-range scalar potential, which suggests that 
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these operators may play a role even in the presence of waves which 
are not unidirectional. 

The operator of eq. ( |CT| ) gives a finite result when applied to 
the solutions of II-III; there is no need to renormalize. However, 
the solutions are not eigenfunctions of this operator. A typical so- 
lution involves n integrations dzi over the locations of the n E^(zi) 
operators contained in the wavefunctional, and the ADM operator 
acts as a "lowering operator" , removing one integration. Hence an 
eigenfunction would have to be an infinite sum over wavefunctionals 
of all possible values of n. It is beyond the scope of this paper to 
investigate the finiteness of the norm of such a sum. 



D The transverse Weyl criterion 

1. Classification of Weyl tensors. 

The Weyl tensor is the part of the Reimann tensor which can be 
non-zero even in empty space, and certain of its components induce 
transverse vibrations when inserted into the equation of geodesic 
deviation [J36|. It is therefore a natural object to work with when 



constructing a criterion for the presence of radiation [R7]]. The con 



struction proceeds in two steps. The first step is a straightforward 
mathematical problem: classify Weyl tensors using their algebraic 
properties. In the second step, one uses physical arguments to de- 
termine the Weyl class(es) most closely associated with radiation. 

To begin with the mathematical problem, there are 10 indepen- 
dent real components of the Weyl tensor, and from these one can 
construct 5 independent complex components which have simple 
duality properties. 

Cabcd = [Cabcd + i(^/2eTXYz)e a bmnC™ | n ] /2] (Dl) 
Cabcd = i(^/2eTXYz)eabmnC^ n . (D2) 

Lower case Roman indices global; upper case Roman 

indices A,B,C,. . . are local Lorentz. e a b mn is the totally antisym- 
metric global tensor, while €txyz is the corresponding local Lorentz 
quantity, the Levi-Civita constant tensor. The duality eigenvalue is 
(5/exxYZ = ±1- There is another Levi-Civita tensor hidden in the 
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^abmn; 

A B C D 

Cabmn = e a e b e c e d €abcd; 

therefore exxYZ and its associated sign convention drop out after 
the conversion to Ashtekar variables and the 3 + 1 splitup. The 
final 3+1 Hamiltonian contains only the phase 5. My convention is 
5 = +1, but in this paper all results are stated in a manner which 
facilitates a conversion to the opposite convention. Of course the 
combinations with simple duality properties also have simple trans- 
formation properties in the local Lorentz frame, which is why one 
chooses to work with C rather than C, when attempting a classifi- 
cation. 

Petrov was the first to classify Weyl tensors by their algebraic 



properties [11]; but for present purposes the equivalent classification 



scheme due to Debever 12, O] is more convenient. A null vector k 



is said to be a principal null vector (Debever vector) of C if 

k [a C b]mn[c k d] k m k n = 0; k a k a = 0. (D3) 

Debever proved that a Weyl tensor can have up to four distinct 
principal null vectors, and he classified Weyl tensors by the number 
of degeneracies among these vectors. If [1111] denotes the Weyl 
tensors which have 4 distinct Debever vectors, [211] the Weyl tensors 
with two vectors degenerate and the rest distinct, etc. , then the five 
classes are [1111], [211], [22], [31], and [4]. (The corresponding five 
Petrov classes are I, II, D, III, and N, respectively.) 

Now suppose k is a Debever vector obtained by solving eq. ( P3|) . 
Make it one leg of a null tetrad k, 1, m, m. Choose the Z axis of a 
local free fall frame so that k and 1 have spatial components along 
±Z, while m and m are transverse. 

- k a l a = m a m a = 1; 
k a k a = l a l a = m a m a 

= k a m a = l a m a = 0. (D4) 

C may be expanded in this basis, and (not surprisingly) one gets 
five possible terms. 

Cabcd = CiV a bV c d 
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where 



+C 2 (V ab M cd + M ab V cd ) 

+C 3 (M ab M cd - U ab V cd - V ab U cd ) 

+C 4 (U ab M cd + M ab U cd ) 

+C 5 U ab U cd , (D5) 



V ab = 2k [a m b ] ; 
U ab = 21[ a m b ] ; 

M ab = 2k [a l b] +2m [a m b] . (D6) 



The five combinations in eq. ( P5|) are the only ones allowed by the 
duality convention 8 — +1. The expansion for the opposite dual- 
ity convention may be obtained from eqs. ( |D5| ) and ( |D6| ) by inter- 
changing m and m in eq. (|D6| ). Eq. ( ]D5| ) is essentially the expan- 



sion given by Szekeres [36], after a relabeling of the basis vectors 
(k, 1, m, m)— *(k, — m, t, t). Since the expansion treats k and 1 quite 
symmetrically, it is valid also for the case that 1, rather than k is 
the principal null vector. 

At this point one turns from the mathematical to the physical: 
which Petrov/Debever class(es), or which term(s) in eq. ( |D5| ), are 
most closely associated with radiation? Consider first which of the 
five tensors in eq. ( [L)5[ ) distorts a cloud of test particles in the manner 
expected for gravitational radiation. Szekeres finds that only the Ci 
and C5 terms produce the transverse displacements in the XY plane 
characteristic of gravitational radiation in the linearized theory. C2 
and C4 produce longitudinal displacements in the XZ or YZ planes. 
C 3 produces a Coulomb (or tidal force) displacement: C 3 distorts a 
sphere of particles into an ellipsoid of revolution with axis along Z. 
These facts suggest that the Ci and C 5 terms signal the presence of 
radiation. 

There is another set of arguments which suggest that the Ci and 
C5 terms are closely associated with the presence of radiation. If the 
Weyl tensor contains only a Ci or C5 term, the tensor is type N. (A 
type N tensor with k as principal null vector contains only a Ci term; 
a type N tensor with 1 as principal null vector contains only a C5 term 
.) Type N is closely associated with radiation. Along characteristic 
curves, when the metric is discontinuous, the discontinuity in the 
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Weyl tensor is type N In the linearized theory, the tensor 

associated with unidirectional gravitational radiation is type N. At 
large distances from bounded sources, the surviving components of 



the Weyl tensor are type N ("peeling theorem" |I3|). 



Although the Ci and C5 terms are closely associated with type N, 
it would be better to call these terms transverse Weyl components, 
rather than type N components, since a tensor which is not type 
N can nevertheless contain Ci or C5 terms. A [22] (type II) field 
contains Ci plus some admixture of longitudinal component, while 
[1111] (type I) contains Ci, C 5 , and C 3 terms. In a theory as non- 
linear as general relativity, one can expect that a collision between 
Ci and C5 transverse waves will produce some C3 (Coulomb) compo- 
nent, and the tensor will be type I rather than type N. In asymptotic 
regions, after the transverse wave has "outrun" its Coulomb com- 
panion, presumably the tensor will revert to type N; but in general 
one should be looking for Ci and C5, rather than type N or any other 
specific Petrov class. One should describe this radiation criterion as 
the transverse Weyl criterion, rather than the type N criterion. 

I now construct operators which project out the Ci and C5 terms. 

d = C abcd l a m b l c m d ; 
C 5 = C abcd k a m b k c m d ; 

C 3 = -C abcd l a m b k c m d . (D7) 

For completeness I have included also the expression for the pure 
Coulomb component C3. The plane wave case has no longitudinal 
components; C2 and C4 vanish identically. 

2. Transition to Ashtekar variables. 

At this point I specialize to the case of plane waves along the 
Z axis. By definition, the plane wave metric has two hypersurface 
orthogonal null vectors which may serve as normals to right- and left- 
moving wavefronts U = (cT-Z)/y / 2 = const, and V = (cT+Z)/a/2 
= const. I identify these normals with k (right-moving) and 1 (left- 
moving), so that a small change in the wave phase will look like 

k a dx a = (-dT + dZ)/V2 = -dU; 

l a dx a = (-dT-dZ)/V2 = -dV. (D8) 
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(Hypersurface orthogonality demands k a dx a oc dU, etc.; the nor- 
malization conditions force the constants of proportionality to be as 
shown in eq. ( |D8| ) ; and the overall phase of k a and l a is fixed by the 
requirement that k° and 1° be positive, i.e. future-pointing.) Lower 
case x denotes a global coordinate; upper case (T, Z, U, V) denotes 
a coordinate in a local Lorentz frame. From the expression for the 
(inverse) tetrads, e A dx b = dX A , k and 1 may be identified with the 
tetrads 



k a = -e aU = +ey; 

l a = e?j. (D9) 



Similarly, 



m a = e a ; 

m a = e a . (D10) 

Evidently the quantities Q are (global scalars and) tensors in a 
Local Lorentz frame. 

The eq. ( [D8|j places quite a strong restriction on the null ba- 
sis, going beyond what is required to maintain the normalization 
eq. (p4|) . The choice eq. ( |D9| ) certainly is not unique. For example, 
k a remains null if it is rescaled by an arbitrary function. (Simul- 
taneously l a must be rescaled by the inverse function, in order to 
maintain the normalization condition — k a l a = 1). Similarly m and 
m may be rescaled. The choice eqs. ( p9| ) and ( |D1Q[ ) facilitates cal- 
culations and leads to highly symmetric formulas for Ci and C5. 
(In this basis, Ci is just C5 with some plus and minus indices inter- 
changed; see eq. ( |1J20[ ) below). For further discussion of the effect 



of choice of basis, see the remarks following eq. ( D20Q . 



Conversion to the Ashtekar language is straightforward. C is 
essentially the (four dimensional) Ashtekar field strength, since the 
C tensor is self-dual, and in empty space the Weyl tensor is the full 
Riemann tensor. 

WFff = e Aa e Bb C abcd . (Dll) 

The four-dimensional field strengths ^F may be replaced by 3+1 
quantities F by using standard formulas. 

< 4 >F™ = -ia5F c M d /2; 

(4) FT = a6 M NsF c s d /2. (D12) 
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M,N,S = space only, cr is a new phase which appears at the 3+1 
reduction step. I choose a = — 1, for reasons explained in Appendix 
A of II. This phase ( unlike 5 ) merely changes the overall sign of 
the Q, and I shall not keep track of the a dependence in the future. 
Useful corollaries of eq. (|D12|) are 

WFff = (i/V2)F±[(5±l)]/2; 

= (VV2)F±[(*=Fl)]/2. (D13) 

If tetrads eqs. (p9|) and ( p!0|) and field strengths eq. (p!3|) are 
inserted into eq. ( p7|) for Ci, the result for 5 = +1 is 

d = i[-F+e° + F+e c z ]e + /2. (D14) 

(For 5 = — 1 replace + by -.) For any metric, typically the Lorentz 
boosts are gauge-fixed by demanding that three of the tetrads van- 
ish: 

e^ = 0, M = space. (D15) 

For the special case of the plane wave metric, the gauge-fixing of 
the XY Gauss constraint and xy spatial diffeomorphism constraints 
imply that four more tetrads vanish. 



■-x — e Y 



0. (D16) 



The tetrad matrix reduces to two 2x2 subblocks which link x,y to 
X,Y (or ±) and z,t to Z,T. Therefore the first term in eq. ( D14j) 
( and only the first term) contains an Ff d term, d = x or y, with 
unacceptable time derivatives of the "coordinate" A~[" . I eliminate 
this term using the classical equations of motion, which are 

(4) F^ d B el = 0, (D17) 

or after 3+1 splitup, and setting B = +, 



= (4)pA- 



cd e A 



- (4) F -+ p c , (4) pU+ c , (4) pV+ c 

— r cd e + ^ r cd e U ^ r cd e V 

= F+e^ 

= F c + d e^ + F+e|. (D18) 

On the second line the ( 4 )F U+ term vanishes because of eq. ( |U 1 3| ) , 
and the ^F^+e^ may be dropped because at the next step the entire 
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term will be contracted with e+. When eq. ( |L) 1 8| ) is inserted into 
eq. ( pi4| ), the result is 



d = iF+e c z e*. (D19) 

This is not quite Ashtekar form, because the triads must be den- 
sitized. Also, the c index can equal z only, because of the gauge 
conditions eqs. ( P15| ) and ( P16|) . The final result is (for C 5 and C 3 
also, since they are calculated similarly) 

Ci = iF z + d E^/( 2 )E ; 
C 5 = iF z - d F//(2)E 

C 3 = F* y /2E|. (D20) 

( 2 *E is the determinant of the 2x2 XY subblock of the tetrad matrix. 
The results for 5 = — 1 are the same, except for overall phases, and 
interchange of + and - everywhere. 

In the case that the wave is unidirectional, the results eq. ( |L)20| ) 
are consistent with the BPR constraints. For example, if the wave 
is right-moving, then the principal vector is k, associated with the 
tensor Ci. From eq. fll4]), A a vanishes, implying that (Ci is finite, 
while) C5 vanishes. 

The Ci of eq. ( p20| ) were calculated in a specific basis; in the lan- 
guage of section III, they are not even kinematic, much less physical. 
In particular the factors of ( 2 >E in eq. ( p20[ ) are basis dependent. 



For example, suppose one shifts from the tetrad basis, eq. (|D9|), to 
a basis in which k a is affinely parameterized. (In the tetrad basis 
one has k b k a; b = Ak a , A 7^ 0.) Then the factor in Ci disap- 
pears, replaced by a factor of E|. One could continue to rescale Ci 
in this manner, until it became density weight unity; then it could 
be sandwiched between holonomies and integrated over z to make it 
kinematic. Presumably one would have to tolerate some degree of 
non-polynomiality in the final result. 

Although it would not be hard to make the transverse opera- 
tor Ci kinematical, it is unikely that Ci could be made physical. 
Gravitational radiation is closely identified with transversality only 
in the the linearized theory. In the full classical theory, scattering 



of two transverse waves produces a C3 Coulomb component [16 . 
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Presumably, then, the commutator of any purely transverse opera- 
tor with the Hamiltonian will not be especially simple, even in the 
classical theory. This is one reason why the main body of the paper 
concentrates on the BPR operators, rather than the Weyl tensor. 
Although it is unlikely that any transverse criterion could be made 
physical, a kinematical criterion for transversality should be both 
feasible and useful. 
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